Relating computational complexity and quantum spectral complexity 
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It is found that the statistical level fluctuations of the AQC 3-SAT problem undergo a transition 
from a poisson (regular) fluctuation form to a form consistent with the predictions of Random 
Matrix Theory. We present data which suggests this transition correlates with the computational 
phase transition in the classical 3-SAT problem. Application to Gaussian Processes and implication 
for experiment is discussed. 
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INTRODUCTION 

The statistical examination of level fluctuation proper- 
ties in the spectra of quantum systems has a rich history 
beginning with the initiation of Random Matrix The- 
ory (RMT) by Wigner [l[ and developed by Dyson 0| 
and Mehta [1] to predict the statistical properties of nu- 
clear spectra and later applied to spectra from atomic 
and condensed matter systems. Such spectra arc charac- 
terized notably by level repulsion and a strong correla- 
tion between energy levels and this earlier work took as 
essential the condition that one is dealing with a system 
containing a large number of degrees of freedom (many- 
particle systems). Later, the domain of validity of RMT 
was enlarged to include the spectra of simple quantum 
systems possessing classical analogs characterized by ir- 
regular (chaotic) motion Q Q . Subsequent research led 
to the correlation between spectral properties predictable 
by RMT and the presence of positive Liapunov expo- 
nents in the corresponding classical dynamics in a large 
number of quantum systems possessing a chaotic clas- 
sical analog. Recently, the presence of level fluctuation 
properties consistent with RMT predictions was found 
in the spectra of Adiabatic Quantum Computing (AQC) 
algorithm for the 3-SAT boolean satisfiability problem, 
an NP-Complete problem Q . 

In this letter, we examine the correlation between spec- 
tral complexity and classical computational algorithmic 
complexity in the AQC 3-SAT, a boolean satisfiability 
problem containing m clauses in n variables, with each 
clause containing 3 boolean variables. The classical com- 
putational complexity of this problem has been studied as 
a function of the clause-to- variable parameterQjll, and 
is characterized by a computational phase transition sep- 
arating computationally 'easy' and 'hard' regions. The 
spectral complexity is measured by the degree to which 
the corresponding adiabatic quantum computation spec- 
trum takes on the spectral properties predicted by Ran- 
dom Matrix Theory. In the following, we present initial 
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results relating the classical and quantum measures of 
complexity for the 3-SAT problem. 



BACKGROUND 

Central to the adiabatic quantum computational 
(AQC) algorithm [| is the adiabatic evolution of the 
ground state energy level, with this state evolving from 
a easily-constructable state to a final state that encodes 
the solution to a particular instance of the 3-SAT prob- 
lem. Until recently, little attention was paid to the 
characteristics of the entire spectrum. However, such a 
characterization revealed the presence of level repul- 
sion and fluctuation properties consistent with RMT in 
portions of the interpolation region for select problem 
parameters. 

This presence of such spectral statistics enlarges fur- 
ther the domain of validity of RMT to include Quan- 
tum Computational systems. It remains to characterize 
the quantum problem via spectral regularity throughout 
problem parameter space. 

Additionally, there has been extensive theoretical and 
experimental work examining systems that obye a Ran- 
dom Matrices type theory containin g an additional adi- 
abatic parameter dependence Such systems, 
called Gaussian Processes, have a number of interesting 
properties including a universal scaling that brings two- 
point correlation functions of eigenvectors or eigenval- 
ues from differing spectra to a common, universal form. 
Thus, the characterization of spectra as adhering to RMT 
would provide the first step in revealing possible proper- 
tics associated with Gaussian processes in AQC. 

This work parallels the correspondence between the 
global change in spectral complexity and classical dynam- 
ical order as measured by the Liapunov dynamical expo- 
nent central to the study of quantal chaos in the 1980s 
and 1990s 0101; extending ideas regarding the Bohr Cor- 
respondence Principle as applied to classically chaotic 
systems. 

This paper is organized as follows: we begin with an ex- 
amination of Random Matrix Theory and spectral mea- 
sures followed by a discussion of the quantum adiabatic 
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computational algorithm. This is followed by the results 
of the statistical analysis of direct diagonalization of the 
AQC 3-SAT spectrum for various problem parameters. 
These results and their generalization to experiment are 
discussed in the context of the asymptotic large N (qubit) 
limit. 



RANDOM MATRIX THEORY 

Random Matrix Theory (RMT) 0] has been useful 
in describing the statistical fluctuations of spectra in 
a variety of quantum systems. Historically, RMT 
has found application describing statistical spectral char- 
acteristics of a variety of complex systems such as 
compound- nuclear resonances [ij, many-particle atomic 
spectra, weakly disordered systems Q , few-body atomic 
systems possessing a chaotic classical analog [ll 15|, and 
even number theory in the distribution of the zeros of the 
Riemann zeta functionlldl 17|[l8|. In addition, recent 



work[6[ has shown the additional applicability of RMT 
in selected instances from the AQC 3-SAT problem 

RMT is a statistical description of complex quantum 
systems in which detailed knowledge about particle in- 
teractions is abandoned in favor of a general ensemble 
description with the only restriction due to physical sym- 
metries. Systems conserving time reversal symmetry arc 
described by the Gaussian orthogonal ensemble (GOE), 
while systems with broken time reversal symmetry arc 
described by the Gaussian unitary ensemble (GUE). 

When examining the global spectral character of a sys- 
tem, the average level density p{E), a model dependent 
property, is normalized to unity. When so normalized, 
the resulting unfolded spectra assumes universal prop- 
erties, and the spectra from a large variety of physical 
systems may be compared, with theoretical predictions 
determined from the theory of Random Matrices. In par- 
ticular, the theoretical spectral nearest-neighbor spacing 
(NNS) fluctuation distribution assumes a universal form 
that may compared with the statistical properties of the 
experimental levels of interest. Whether the specific dis- 
tribution assumes universal properties depends on the 
presence of quantal symmetries in the system, and, in 
the case of quantal systems with classical analogs, RMT 
spectral predictions correlate with the presence of classi- 
cal nonlinearities as distinguished, for example, by posi- 
tive Liapunov exponents. 

The NNS distribution typical for complex quantum 
and disordered systems is given by the Wigncr distribu- 
tion d, as predicted by RMT. Such distributions are 
typical for complex systems without quantal symmetries, 
which results in correlated, irregular energy spectra that 
exhibit level repulsion. Conversely, Hamiltonians corre- 
sponding to classical systems subject to symmetries and 
conservation laws typically display uncorrelated energy 
spectra and nearest-neighbor level spacing distribution 



p(x) that is Poisson-like {p{x) ~ exp{—x)). 

The spectral regularity may be quantified by the Brody 
parameter q [l^ . occurring in a one-parameter distribu- 
tion that interpolates between a regular Poisson spec- 
trum {q = 0) and an irregular Wigner distribution 
(q = 1) [l9| . A renormalized spectrum having Brody 
parameter q is characterized by the following NNS prob- 
ability distribution for level spacing 6, for the GOE: 

P,{S) = {1 + q)f3S'' exp(-/3<5i+«) , /3 = T (^^^ . 

(1) 

The AQC formulation is defined in terms of an adi- 
abatic parameter that interpolates between an initial 
and final Hamiltonian. In the theoretical limit of in- 
terpolation between two independent random matrices, 
parametric random matrix theory, or Gaussian processes 
(CP) is invoked. Thus, one parameter disordered spectra 
is compared with the theoretical predictions of Gaussian 
processes (GP) whose two point correlation function as- 
sume a universal form when the parameter x is rescaled 

by the average spectral quantity ei, a; — > \/ ( ^)'^ x as 
shown in the context of parametric RMT (GP) and 
derived also in the context of supersymmetry models 
Such spectra are characterized by an abundance 
of avoided level crossings, such as the spectra shown for 
the present case of the AQC 3-SAT problem shown in 
Fig. [H Additional universal results have been obtained 
for multiparameter GP that relates to the occurrence of 
Berry phase in such systems and the (constant) universal 
density of Berry points in such systems when rescaled as 
above. 

The GP in one parameter can be defined in a manner 
similar to RMT through its first two moments 



H,j{x) 



H,j{x)Hu{x') 
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where g\l^P 



J/3=2) 



Sik^ji + 5ii5jk, gtj^ki' = '^SiiSjk, u) IS re- 
lated to the Hamiltonian dimension and mean level spac- 
ing, and f(x) is normalized such that /(O) = 1. (3=1 
is the GOE case, consisting of H{x) that are real sym- 
metric matrices, and /3 = 2 is the GUE case consisting of 
H{x) that are complex Hermitian matrices. 

A natural extension of RMT, Gaussian processes have 
been examined in the context of quantum dissipation 

mill [13, with a major component of the dissipation re- 
sulting from Landau-Zener nonadiabatic level transitions 
24 1 at avoided level crossings [2^. This aspect is impor- 
tant to the adiabatic quantum computation algorithm, as 
algorithmic success requires the evolution of the system 
while remaining in the ground state, without such tran- 
sitions. As a result, interpolation through irregular spec- 
tral regions have a higher probability of transition and 
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the interpolation must be slowed in this region when em- 
ploying this approach to prevent level transitions. From 
a computational standpoint, algorithmic cost increases 
with spectral complexity characterized by the onset of 
random matrix (Wigner) statistics. For an AQC system 
possessing uniform average level spacing and obeying the 
conditions of a GP, it can be shown through the theory 
of GP that a quantum exponential speedup over classi- 
cal algorithms is unlikely [6|; however, such a result for 
the more realistic physical condition of a level density 
varying with energy is currently lacking. 

THE QUANTUM ADIABATIC COMPUTATION 
ALGORITHM 

The adiabatic algorithm allows for a quantum adi- 
abatic process to perform a computation of an NP- 
Completc problem, the 3-SAT boolean satisfiability prob- 
lem. The time dependent Hamiltonian for such a process 
with an interpolation time T takes the form 

H{^)^{l-^)H, + ^H,, (3) 

where the Hamiltonian Hjj has the property that the 
ground state is easily constructed, and the Hamiltonian 
Hp has the property that the ground state encodes the 
solution to the particular instance of the 3-SAT problem 
under investigation. 

Starting with the ground state eigenstate of i/b, the ap- 
proach allows for the calculation of a nontrivial computa- 
tion problem through adiabatic interpolation of this state 
from the initial to final Hamiltonian. Upon completion of 
the ground state interpolation, the state encodes the so- 
lution to the 3-SAT problem, which can be read through 
projective measurement of the state. This approach has 
been shown to be equivalent to 'traditional' circuit mod- 
els of quantum computation [2^, in which quantum algo- 
rithmic speedup over classical approaches for a number 
of algorithms is well established. 

The traditional measure of algorithmic complexity of 
the 3-SAT problem using classical computational mea- 
sures involving either the Davis-Putnam 0] or the GSAT 
algorithm Q , shows a marked transition at a given clause- 
to-variable ratio, / ~ 4.2, indicating a computational 
phase transition from easy to hard solubility [27[ in the 
large n limit. For problem sizes having n < 20, problem 
sizes typical in numerical simulation, the transition phe- 
nomenon still exists, though the transition point is less 
well defined and somewhat displaced from its asymptotic 
value. 

RESULTS 

The statistical spectral fluctuation analysis involved 
determining the NNS distribution of the Hamiltonian 
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FIG. 1: Spectrum of irregular unfolded energies of the in- 
stantaneous 3-SAT quantum computational Hamiltonian for 
N=8 m=16 clauses in the A'^ ~ [40,60] region. The predom- 
inance of avoided level crossings is characteristic of complex 
(irregular) spectra. 

H{t). Eq. ([3]), at 100 points in the interpolation of a 
3-SAT Hamiltonian with n = 8 variables having matrix 
dimension N = 256 for several values of the clausc-to- 
variable ratio /. In light of the above discussion relat- 
ing computational cost and spectral complexity, we are 
interested in determining whether an irregular spectral 
region occurs in the interpolation and for each problem 
instance, the maximal Brody parameter is determined, 
with the results averaged over 200 problem instances for 
each /. 

The results reveal a systematic change in the spec- 
tral regularity (Brody measure q) of the instantaneous 
AQC Hamiltonian H(t) during the course of the adiabatic 
quantum algorithm. A typical unfolded (renormalizcd) 
spectrum associated with a particular problem instance 
is shown for a particular problem instance having N = 8, 
/ = 2, shown in Fig. [TJ In the initial phase of the inter- 
polation for hard-to-solve problem instances, the statisti- 
cal NNS fluctuations conform to a regular, Poisson-type 
distribution. This corresponds to eigenvalue degenera- 
cies inherent in the initial hamiltonian, H[0). Later in 
the interpolation, the fluctuation probability distribution 
transitions to an irregular, Wigner-type distribution. It 
was found that such irregular spectra only occur for com- 
putational problem instances having (/ > 1). 

Fig. [2]summarizes these quantum simulation results as 
a function of the classical algorithmic complexity param- 
eter /, plotting spectral complexity (Brody measure) ver- 
sus / (clause-to- variable ratio). For easy-to-solve compu- 
tational problem instances (/ small), the spectral com- 
plexity was found to be zero. For hard-to-solve instances 
of the 3-SAT problem the spectral complexity is nonzero. 
The correlation between the "classical" computational 
complexity as parameterized by / and the quantum spec- 
tral complexity, parameterized by the Brody parameter 
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FIG. 2: Spectral Complexity, maximum Brody parameter, as 
a function of the classical computational complexity, clause- 
to-variable ratio, of the 3-SAT problem for N=8. 



q is shown in Fig. [21 Wc note that if the computational 
complexity is used as a proxy for computational cost, 
the quantum system appears to go through an orderly- 
complex transition in the region 1 ^ / ^ 4 and dimin- 
ishes somewhat for f ^ 5. The overall behavior exhibits 
the basic features of the classical algorithmic cost pro- 
file for the 3-SAT, which, for n = 8, is characterized by 
a computational phase transition that is smoothed com- 
pared to the asymptotic form and displaced from the 
asymptotic (large n) phase transition point of 4.2. 

Additional studies for n closer to the asymptotic re- 
gion {n > 20) would be required to examine the asymp- 
totic quantum transition and comparing the classical and 
quantum cost profiles in detail. While the current results 
arc for the n ~ 8 non- asymptotic regime, the change in 
global spectral statistics correlated with computational 
complexity as measured by clause-to-variable ratio was 
found to be significant. The salient feature is the tran- 
sition to from regular to irregular spectra dependent on 
algorithmic complexity that correlates with classical re- 
sults. 



DISCUSSION 

Recent work [§| showed that an NP-Complete problem 
can be converted into an adiabatic quantum mechanical 
process, which marked a new approach to quantum com- 
putation that was later shown equivalent to the circuit 
model of quantum computation[26[. 

The present results enlarge the regime of validity of 
RMT to this additional class of AQC problems. The cen- 
tral question we ask in this article is whether the spectral 
complexity is correlated with computational algorithmic 
complexity. 

The results from Fig. [2] illustrate that when the entire 
spectrum from the AQC problem instances are examined. 



the spectrum is shown to be orderly for small clause-to- 
variable ratio, /, becoming more complex for larger /. 
This leads us to the following result: when formulated 
in terms of an adiabatic quantum computation, all com- 
putationally hard NP-complete complexity class 3-SAT 
problems traverse a region characterized by an irregular 
spectrum. 

How general is this result? Our results verify the above 
conclusion for n = 8 in regions that have been character- 
ized as algorithmically complex (/ ~ 4.2), and we sug- 
gest that this result likely holds for larger n. Spectra 
exhibiting orderly spectral statistics (g ~ 0) are typically 
associated with quantum symmetries, and it seems un- 
likely that a symmetry not manifesting itself for n = 8 
will manifest itself for larger n. 

The above conjecture can be tested through a related 
experimental hypothesis related to the spectral statis- 
tics of the AQC problem: there exists no classically hard 
instance of an NP-complete complexity class problem 
which, when formulated in terms of an adiabatic quan- 
tum computation, results in a simple spectrum (Poisson 
level spacing distribution) at all interpolation points in 
the computation. 

Further, since quantum adiabatic computing has been 
shown equivalent to the circuit model of quantum com- 
puting (based on quantum bits) [3], the above predic- 
tions would apply equally to all quantum computation- 
ally hard circuit model problems as well. 

The experimental determination can be made by ex- 
amining the spectrum at intermediate points in the inter- 
polation for problem instances having fixed n and vary- 
ing clause-to-variable ratio. This hypothesis would be 
disproved if a system were found that has the properties 
of (1) being NP-Complete, (2) computationally complex 
(hard-to-solve instance) determined by classical measures 
such as clause-to-variable ratio, and (3) the property of 
having a regular spectrum (Brody parameter g = 0) at all 
interpolation points. Such a system would correspond to 
an NP complete, hard-to-solve problem instance possess- 
ing a regular spectrum, which would require no inter- 
polative slow down due to irregular (complex) spectral 
effects. 

Finally, we have identified hard-to-solve instances of 
NP-Complete quantum adiabatic computation processes 
as having parameter regions obeying RMT and thus be- 
longing to the class of Gaussian processes (CP). Once this 
identification is made, a variety of universal results fol- 
low relating primarily to two-point correlation functions 
of wavcfunctions and eigenvalues that have been tested in 
quantum systems possessing a classical analo g an d con- 
taining an external parameter dependence (lot .[iol . If the 
above conjecture happens to be true, then spectral com- 
plexity would be correlated with computational algorith- 
mic complexity and results consistent with CP should be 
found asymptotically in AQC processes having a problem 
parameter / consistent with significan computational al- 
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gorithmic complexity. 
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